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ABSTRACT 

Snap  roll  is  a  phenomenon  which  has  plagued  the  sub- 
marine designer  since  the  advent  of  the  high  speed  submarine 
A  description  of  snap  roll  and  its  possible  causes  is  pre- 
sented in  this  thesis.   A  computer  program  is  developed 
which  simulates  a  submarine  in  surge,  sway,  yaw  and  roll. 
It  is  verified  that  snap  roll  is  very  sensitive  to  the 
metacentric  height.   Additionally,  it  is  seen  that  reducing 
the  sail  size  does  not  have  as  large  an  effect  on  reducing 
the  roll  angle  as  anticipated.   Rudder  sequencing  and  speed 
reduction  are  two  other  methods  examined.   The  results  of 
the  four  investigations  are  presented  in  graphic  form.   The 
conclusions  include  a  justification  for  the  use  of  the 
simulation  computer  model  as  a  design  tool. 
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NOMENCLATURE 

IYMBOL  DEFINITION 

CB  Center  of  buoyancy. 

CG  Center  of  gravity. 

D  Drag  force. 

I  Moment  of  inertia  about  the  i  axis 


I.,  Product  of  inertia  about  jk  axis. 

K  Hydrodynamic  moment  about  x  axis 

(rolling  moment) . 


K.  Non-dimensional  coefficient  used  in 

1  representing  K  as  a  function  of  i  . 

0  Overall  ship  length. 

L  Lift  force. 

m  Ship's  mass. 


N  Hydrodynamic  moment  component  about 

z  axis  (yawing  moment). 


N  .  Non-dimensional  coefficient  used  in 


representing  N  as  a  function  of  i 
Angular  velocity  about  the  x  axis 


p  Angular  acceleration  about  the  x  axis 

q  Angular  velocity  about  the  y  axis. 

q  Angular  acceleration  about  the  y  axis 

r  Angular  velocity  about  the  z  axis. 

r  Angular  acceleration  about  the  z  axis 


U  Velocity  of  the  origin  of  the  body  axes 

relative  to  the  fluid. 


u  U  velocity  component  in  the  x  direction 


u  Acceleration  component  in  the  x 

direction. 


u  Initial  steady  value  of  U  in  the  x 

direction. 

v  U  velocity  component  in  the  y  direction 


Acceleration  component  in  the  y 
direction. 


w  U  velocity  component  in  the  z  direction 


w  Acceleration  component  in  the  z 

direction. 


xR  The  x  coordinate  of  the  center  of  buoyancy 

xn  The  x  coordinate  of  the  center  of  gravity. 

G 
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x  The  x  coordinate  of  the  center  of 


s 


pressure  of  the  sail 


Hydrodynamic  force  along  x  axis 
(axial  force) . 


X '.  Non-dimensional  coefficient  used  in 

1  representing  X  as  a  function  of  i. 


B 


s 


The  y  coordinate  of  the  center  of 
buoyancy. 


y_  The  y  coordinate  of  the  center  of 

gravity. 


y  The  y  coordinate  of  the  center  of 


pressure  of  the  sail. 


Hydrodynamic  force  along  y  axis 
(lateral  force) . 


Y*.  Non-dimensional  coefficient  used  in 

1  representing  Y  as  a  function  of  i . 


zR  The  z  coordinate  of  the  center  of 

buoyancy. 


zr  The  z  coordinate  of  the  center  of 

gravity. 


z  The  z  coordinate  of  the  center  of 

pressure  of  the  sail. 


A  Angle  of  drift. 

oh  Bowplane  or  sailplane  deflection. 
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O  Rudder   deflection 


r 


5  Sternplane  deflection 


s 


/  Ratio  u  /u. 
1  o/ 

"  Angle  of  pitch 

•  Angle  of  yaw. 

I  Angle  of  roll. 


a-  ,  b.  ,  c  Constants  used  to  represent  the 

1   1  propeller  thrust  in  the  axial 

equation. 
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CHAPTER  I  -  INTRODUCTION 

1.1   BACKGROUND 

Since  the  advent  of  the  high  speed  submarine  in  the 
late  fifties,  the  submarine  designer  has  been  confronted 
with  new  design  problems  in  the  area  of  ship  control.   One 
of  these  problems,  which  is  still  with  us  two  decades  later, 
is  snap  roll.   Part  of  the  reason  for  the  perpetuation  of 
this  problem  has  been  a  general  lack  of  adequate  data  and 
information  concerning  this  phenomenon.   Consequently, 
appropriate  design  criteria  could  not  be  developed  and  the 
problem  of  snap  roll  has  continued.  Naval  architects  and 
ship  designers  need  adequate  alternatives  and  trade-offs  to 
provide  good  designs.   This  is  the  motivation  for  this 
thesis. 

A  computer  model  will  be  developed  in  Chapter  II  which 
simulates  a  submarine  in  roll,  yaw,  surge  and  sway.   This 
model  can  be  used  by  naval  architects  as  a  design  tool  in 
performing  trade-off  studies  or  as  a  method  of  establishing 
design  criteria  for  a  particular  design.   The  alternatives 
for  reducing  snap  roll,  which  are  investigated  in  Chapter  III 
using  this  simulation  model,  were  based  on  the  immediate 
needs  of  the  design  community  and  the  time  and  monetary 
constraints  normally  placed  on  this  type  of  effort.   It  was 
felt  that  these  alternatives  should  not  increase  the  com- 
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plexity  of  present  designs  or  alter  present  design 
practices,  in  order  that  they  may  be  put  to  immediate  use. 
Chapter  IV  presents  the  conclusions  and  recommendations  of 
this  thesis.   In  the  appendix  will  be  found  a  listing  of 
the  computer  model,  a  flow  chart  and  a  sample  output. 

1.2  WHAT  IS  SNAP  ROLL? 

An  excellent  description  of  snap  roll  was  presented  by 
Griffin,  et  al  (l)1: 


"At  20  sec  after  full  rudder,  the  initial 
transients  have  died  out,  and  the  lateral 
states  begin  to  follow  the  trim  values 
which  correspond  to  the  decreasing  forward 
speed.   Eventually,  the  forward  velocity 
reaches  its  trim  value,  and  r,  q,  v,  w 
and  *f  are  fixed  at  their  respective  trim 
values.   For  r,  q,  v,  and  w,  the  difference 
between  the  plateau  value  and  the  later  trim 
value  is  not  large;  however,  for  roll,  it 
can  differ  by  more  than  13  deg.   This  is 
snap  roll. " 


Simply  stated,  snap  roll  is  the  maximum  roll  amplitude 
that  results  from  a  turning  maneuver.   The  large  roll  angle, 


Numbers  in  parentheses  indicate  references  listed 
at  the  end  of  this  thesis 
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which  develops  early  in  the  turn,  constitutes  the  "snap." 
The  large  peak  value  of  the  roll  angle  is  the  submarine 
tending  towards  its  steady  state  turning  value  for  the 
initially  high  forward  speed.   This  usually  occurs  about 
20  seconds  into  the  turn.   As  the  submarine  enters  the 
turn  and  during  the  initial  phase  of  the  turn,  there  is  a 
large  side  slip  velocity.   This  side  slip  velocity  probably 
accounts  for  a  major  portion  of  the  maximum  roll  angle. 
The  subsequent  reduction  of  the  large  roll  angle  is  pri- 
marily due  to  speed  loss.   The  high  drag  force  developed 
during  the  turn  maneuver  dramatically  reduces  the  forward 
velocity  component.   The  steady  state  forward  speed  can  be 
expected  to  be  1/4  to  1/3  of  the  initial  speed. 

There  is  still  conjecture  amongst  hydrodynamicists  as 
to  the  theory  and  physics  involved  in  explaining  this 
phenomenon.   It  is  generally  agreed  that  snap  roll  pri- 
marily depends  on  the  turning  rate,  the  amount  of  rudder 
deflection,  the  fairwater  (sail)  size  and  location,  and  the 
ship  speed.   A  large  turning  rate  and/or  rudder  deflection 
produces  a  large  lateral  velocity.   Consequently,  this  large 
v  acting  on  the  sail  produces  a  large  roll  moment.   Addi- 
tionally, it  has  been  proposed  that  the  starting  lift 
exceeds  the  steady  lift  due  to  the  dominating  effect  of  the 
trailing  vortex  sheet  downwash  in  the  steady  case  (2). 
This  may  be  due  to  the  starting  vortex  formed  and  shed  by 
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the  sail.   However,  there  is  no  substantiating  evidence  for 
this. 


m- 


CHAPTER  II  -  THE  SIMULATION  MODEL  AND 
COMPUTER  PROGRAM 


II. 1   THE  FOUR  DEGREE  OF  FREEDOM  MODEL 

There  exists  at  least  three  computer  programs  which 
simulate  a  submarine  in  a  maneuver.   However,  none  of  these 
programs  are  available  for  general  use,  primarily  due  to 
proprietary  considerations.   Therefore,  in  order  to  inves- 
tigate snap  roll  in  submarines,  it  was  necessary  to  develop 
another  simulation  model.   Unlike  the  other  models,  however, 
this  model  involves  the  use  of  four  degrees  of  freedom  vice 
six  degrees  of  freedom  equations.   The  four  coupled  equations 
considered  here  are  the  X,  Y,  N  and  K  equations.   The  justi- 
fication for  the  use  of  only  four  equations  of  motion  is 
that  in  a  turn,  depth  excursions  are  small  and  the  Z  and  M 
equations  may  be  decoupled.   The  depth  change  which  occurs 
during  a  turn,  called  "squat",  normally  can  be  easily 
remedied  with  a  small  deflection  of  the  diving  planes  and/or 
sternplanes.   These  small  excursions  have  a  negligble  effect 
on  snap  roll.   This  is  evidenced  by  the  fact  that  a  four 
degree  of  freedom  simulation  program  provides  valid  results 
for  a  turning  maneuver.   This  will  be  seen  from  the  results 
presented  later  in  this  chapter  and  in  Chapter  III.   The 
advantages  afforded  by  a  four  degree  of  freedom  program, 
in  addition  to  providing  valid  results,  are: 
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•  Reduced  computation  time. 

•  The  estimation  or  measurement  of  fewer 
hydrodynamic  coefficients. 

Of  course,  such  a  model,  necessarily  restricts  one's  con- 
sideration to  only  maneuvers  involving  surge,  sway,  roll 
and  yaw. 

II. 2   THE  COORDINATE  SYSTEM  AND  EQUATIONS 

The  coordinate  system  chosen  for  use  in  this  model  was 
based  on  reference  (3).   This  reference  gives  the  standard 
equations  of  motion  for  submarines  used  by  the  U.S.  Navy. 
The  coordinate  system  is  illustrated  in  figure  2-1.   It 
should  be  noted  that  the  origin  is  centered  at  the  center  of 
gravity.   A  better  choice  may  have  been  to  use  the  center  of 
buoyancy;  however,  to  avoid  a  rewriting  of  already  valid 
equations,  the  reference  (3)  system  was  used.   The  center  of 
buoyancy  of  the  basic  hull  for  a  submarine  design  is 
necessarily  the  geometric  center  of  the  hull.   Therefore,  a 
transverse  velocity  (or  acceleration)  would  not  produce  a 
roll  moment,  caused  by  the  basic  hull  form  and  this  would 
reduce  the  number  of  necessary  calculations. 

The  equations  used  in  the  development  of  this  simu- 
lation model  were  taken  directly  from  reference  (3)  and  are 
listed  in  Appendix  B.   There  are  no  second  order  or  higher 
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order  acceleration  terms  involved  in  these  equations.   This 
is  due  to  the  lack  of  any  significant  interaction  of  the 
viscous  and  inertia  forces.   This  allows  the  use  of  potential 
theory  which  gives  adequate  values  for  the  hydrodynamic 
forces  (4).   Acceleration  and  velocity  interaction  terms 
are  also  assumed  negligibly  small  for  the  same  reason  stated 
above.   It  should  be  noted  that  the  equations  remain  non- 
linear even  without  these  second  order  acceleration  terms, 
due  to  the  higher  order  velocity  terms. 

The  equations  are  applied  to  a  submarine  in  a  turning 
maneuver  under  the  following  assumptions  and  restrictions: 

1.  The  submarine  is  neutrally  buoyant,  i.e. 
W  -  B  =  0. 

2.  The  submarine  is  at  least  three  hull 
diameters  below  the  free  surface,  i.e. 
fully  submerged. 

3.  The  z  component  of  the  velocity  is  zero. 

4.  There  is  no  initial  list,  i.e.  f*0  -  0. 

5.  There  is  zero  trim,  i.e.  9  =  q  =  q  =  0. 

6.  The  origin  is  at  the  center  of  gravity, 
i.e.  xQ  =  yG  =  zG  =  0. 

7.  Only  the  moments  of  inertia  along  the 
principal  axes  are  non-zero,  i.e. 

Ixy  =  Iyz  =  Ixz  =  °"   This  is  Partly  a 
consequence  of  #6. 
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8.  The  coefficients  and  equations  used  are 
in  non-dimensional  form. 

9.  The  control  surface  deflections  are 
initially  zero,  with  only  the  rudder 
angle  being  non-zero  for  t>0. 


II. 3   THE  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

The  simplest  method  of  solving  for  the  roll  angle 
induced  by  a  turning  maneuver  is  to  solve  for  u(t)  from 
the  non-linear  X  equation;  then,  to  use  this  u  in  solving 
the  linear  Y,  N  and  K  equations.   This  will  provide  com- 
paratively accurate  values  for  the  roll  angle  for  relatively 
small  deflections  and  for  small  forward  speed  losses.   How- 
ever, we  are  interested  in  high  speed  and  tight  turns, 
where  snap  roll  is  most  prominent.   The  vehicle  motions  in 
this  type  of  turning  maneuver  are  outside  the  valid  linear 
range.   Therefore,  the  non-linear  equations  of  motion  are 
used. 

The  general  method  used  for  solving  the  four  non- 
linear equations  is  an  iterative  technique.   The  specific 
technique  employed  in  this  thesis  was  adapted  from  the 
method  used  in  reference  (5)-   The  four  equations  are  each 
written  in  terms  of  the  accelerations  u,  v,  r  and  p.   They 
take  the  form: 
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X:   (m'-  X^)  u  -  yQr  =  f1  (u,  v,  r,  p,  £r) 


Y:   (m'-  Yp  v  -  (m'zG  +  Y^  1)  p  + 

(m'xG  -  Y^  1)  =  f2  (u,  v,  r,  p,  Sr) 


n:    (i;  -  Np  r  -  (i;x  +  n;)  p  + 

(m'xG/l2  -  Nt/1)  v  -  (m'yG/l2)  u  = 
f~  (u,  v,  r,  p,  Sr) 

K:   (i;  -  K!)  p  -  (i;z  4  K!)  r  - 

(mfzG/l2  +  Kyi)  v  =  f 4  (u,  v,  r,  p,  Sr) 


The  prime  indicates  a  non-dimensional  value. 

Expressions  for  the  accelerations  u,  v,  r  and 
p  are  derived  by  employing  matrix  methods.   Based  on  the 
initial  conditions  u  ,  v  ,  etc.,  starting  values  for  u,  v, 
r  and  p  are  obtained  for  t  =  0.   Values  for  u,  v,  r  and  p 
are  calculated  from  these  starting  values  employing  the 
Taylor  expansions, 

u(t  +  At)  a  u(t)  +  At-u(t) 

v(t  +  At)  ?   v(t)  +  ^t-v(t) 

r(t  +  At)  -  r(t)  +  At-f(t) 

fit  +At)  -?(t)  +  &t-p(t)  +  ((&t)2).p(t). 

These  new  velocity  values  are  then  used  to  update  the  pre- 
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viously  calculated  accelerations.   This  indicates  the 
iterative  nature  of  this  integration  technique.   The  value 
of  S  and  t  are  updated  after  each  iteration.   The  equations 
are  integrated  over  the  time  span  specified.   The  size  of  At 
is  arbitrary;  however,  the  smaller  the  At,  the  greater  the 
accuracy  obtained  from  this  technique.   The  accuracy  desired 
should  be  tempered  by  the  fact  that  small  At's  require 
rather  large  amounts  of  computer  time.   It  was  found  that 
At's  on  the  order  of  one  half  to  one  second  provided  very 
satisfactory  results,  while  not  using  excessive  amounts  of 
computer  time. 

II. 4   COMPUTER  SIMULATION  MODEL  VALIDITY 

A  large  number  of  computer  runs  were  executed  for 
various  submarines,  under  various  conditions.   It  was  in- 
tended to  prove  the  validity  of  the  developed  computerized 
simulation  model  using  hydrodynamic  coefficients  taken  from 
towing  tank  model  tests.   The  results  from  one  of  these 
computer  runs  is  presented  in  figure  2-2.   As  can  be  seen, 
in  general,  the  model  accurately  predicts  the  submarine's 
response  in  roll;  however,  the  predicted  snap  roll  angle  is 
too  large.   This  can  be  attributed  to  a  number  of  causes. 
One  is  the  possibility  of  scale  effects,  since  the  co- 
efficients used  in  the  program  are  not  the  full  scale 
hydrodynamic  coefficients.   Secondly,  the  metacentric  height 
used  was  that  specified  in  publications.   The  actual  GM 
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during  the  full  scale  trials  could  vary  slightly  from  this 
value  and,  as  will  be  seen  in  Chapter  III,  the  maximum  roll 
angle  is  very  sensitive  to  small  changes  in  GM.   However,  for 
other  submarines  and  sequences,  the  predicted  maximum  roll 
angles  were  less  than  the  full  scale  values  and  therefore, 
this  somewhat  indicates  that  the  difficulty  in  predicting 
the  snap  roll  angle  does  not  lie  entirely  in  the  model  or 
technique  used.   It  will  also  be  noted  that  the  model 
results  increasingly  deviate  from  the  full-scale  trials  as 
t  increases.   This  can  be  attributed  to  two  causes.   First, 
the  full-scale  trials  are  performed  with  a  spiraling 
maneuver.   In  order  to  save  computer  time  and  since  we  v/ere 
only  interested  in  the  maximum  roll  angle,  the  computer 
simulation  runs  were  only  performed  for  180   turns.   There- 
fore, the  rudder  angle  was  taken  off  earlier  for  the  model 
and  the  roll  angle  decreased  more  rapidly  for  large  t. 
Secondly,  at  slower  speeds  (large  t),  the  cross  coupling 
effect  of  pitch  and  heave  on  roll  is  more  pronounced  than  at 
higher  speeds.   Since  we  are  using  only  four  degree  of 
freedom  equations,  the  effect  of  pitch  and  heave  at  low 
speeds  is  lost. 
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CHAPTER  III  -  METHODS  FOR  REDUCING  SNAPROLL 

111.1  GENERAL 

There  are  a  number  of  techniques  and  methods  which  have 
been  proposed  to  reduce  snap  roll  in  submarines  (1).   Four  of 
these  alternatives  will  be  investigated  here.   As  was  in- 
dicated in  the  introduction,  the  choice  of  these  alter- 
natives was  based  on  the  needs  of  the  design  community  and 
time  constraints.   The  four  methods  chosen  can  be  categorized 
under  one  of  two  headings.   The  first  concerns  changes  in  the 
naval  architectural  characteristics  of  the  design.   Increasing 
GM  and  reducing  sail  size  are  classed  under  this  category. 
The  second  category  involves  those  alternatives  which  make 
use  of  an  automatic  ship  control  system.   Rudder  sequencing 
and  speed  reduction  fall  under  this  classification.   Each 
method  of  snap  roll  reduction  will  be  investigated  indi- 
vidually in  the  following  sections. 

111. 2  METACENTRIC  HEIGHT 

The  metacentric  height  is  directly  related  to  the 
righting  moment,  K«,  by  the  relation, 

K«p=  -W-GM 
where  W  is  the  submerged  displacement.   This  term  is  some- 
what disguised  in  the  governing  equations  in  Appendix  A. 
In  the  K-equation, 

24 


K<p  =  zR  •  B 
where  B  is  the  buoyancy,  and  equal  to  the  submerged  dis- 
placement at  neutral  buoyancy,  and  zR  is  the  metacentric 
height.   It  can  be  shown  that,  as  a  first  approximation,  the 
maximum  roll  angle  is  inversely  proportional  to  the  meta- 
centric height  (1).   Therefore,  it  would  be  advantageous  to 
increase  the  metacentric  height  in  order  to  reduce  the 
effect  of  snap  roll.   However,  this  could  prove  to  be  a 
difficult  design  task,  since  it  would  involve  adding  weight 
or  rearranging  weights  in  an  already  weight-limited  design. 

The  metacentric  height  for  several  submarines  was 
increased  by  0.25,  0.50,  0.75,  1-0  and  1.5  feet.   A  non- 
dimensional  plot  of  the  results  are  shown  in  figure  3-1  • 

It  should  be  noted  that  the  amount  of  rudder  used  and  u 

o 

were  different  for  each  vehicle.   This  indicates  that  the 
percent  decrease  in  the  maximum  roll  angle  is  essentially 
independent  of  u  and  the  maximum  ordered  rudder  angle  for 
a  given  GM  change.   The  curve  is  asymptotic  for  greater  than 
100%  changes  in  metacentric  height.   Therefore,  increases  in 
the  metacentric  height  of  greater  than  about  100%  will  yield 
no  greater  than  approximately  an  82%  decrease  in  the  maximum 
roll  angle. 

As  an  example,  a  30%  increase  in  GM  yields  a  33%  de- 
crease in  f*   .   For  a  typical  submarine,  this  translates 
max         J 

into  a  10°  decrease  in  the  maximum  roll  angle  for  an  approx- 
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imately  k   inch  increase  in  the  metacentric  height.   It  can  be 
seen  that  large  decreases  in  the  snap  roll  angle  can  be 
realized  for  relatively  small  changes  in  GM.   This  is 
evident  in  figure  3-2.   This  is  a  plot  for  a  typical  sub- 
marine which  shows  the  decrease  in  the  maximum  roll  angle 
for  increasing  GM.   An  average  5-  75 °  decrease  in  *P    is 
gained  for  each  3  inch  increase  in  GM.   A  point  of  minor 
interest  shown  in  this  plot  is  that  the  maximum  roll  angle 
occurs  earlier  in  the  turn  as  GM  increases. 

III. 3   RUDDER  SEQUENCING 

Rudder  sequencing  means  to  defer  the  full  rudder  angle 
ordered,  until  some  predetermined  speed  loss  has  been  reached 
It  is  expected  that  by  limiting  the  initial  rudder  angle, 
the  maximum  roll  angle  will  be  reduced.   This  is  based  on 
the  fact  that  the  roll  angle  is  directly  proportional  to  the 
rudder  deflection.   Rudder  sequencing  is  envisioned  as  part 
of  an  automatic  ship  control  system;  however,  the  method  of 
rudder  sequencing  is  simple  enough  to  be  employed  as  a 
manual  process.   One  disadvantage  of  this  method  of  snap  roll 
reduction  is  that  we  might  expect  the  transfer  to  increase. 
However,  such  an  increase  is  normally  no  more  than  one  half 
to  one  ship's  length.   The  procedure  employed  in  this  thesis 
to  investigate  this  roll  reduction  method  involves  reducing 
the  maximum  ordered  rudder  angle  by  5»  10 »  15  and  20  degrees. 
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The  predetermined  speed  loss  criteria  is  to  maintain  the 
reduced  rudder  deflection  until  the  speed  in  the  turn  has 
been  reduced  to  one  third  that  of  the  initial  speed  into  the 
turn.   This  is  depicted  in  figure  3-3-   The  selection  of 
this  speed  criteria  was  based  on  the  fact  that  the  steady 
speed  in  a  fully  developed  turn  is  approximately  one  third 
of  the  initial  speed.   The  rudder  deflection  rate  was  held 
constant  for  all  submarines  and  all  sequences  at  U   deg/sec. 

The  results  are  depicted  in  figure  3-^-   As  can  be  seen 
from  the  plot,  this  method  of  snap  roll  reduction  is  effec- 
tive.  However,  relatively  large  rudder  angle  reductions  are 

needed  to  produce  small  changes  in  T       .   For  example,  a  10$ 
^  to       max  ^   ' 

reduction  in  the  maximum  roll  angle  requires  a  23f;  reduction 

in  the  initial  rudder  angle.   The  use  of  rudder  sequencing 

for  a  typical  submarine  is  shown  in  figure  "3-5.   It  can  be 

seen  that  a  greater  decrease  in  T    per  degree  of  rudder  is 

max 

gained  for  the  20  rudder  sequence  than  for  the  30°  rudder 
sequence.   It  should  be  noted  that  some  of  the  figure  3-^ 
data  points  were  based  on  the  same  submarine   at  different 
initial  speeds,  i.e.  2^  and  10  knots.   Therefore,  rudder 
sequencing  seems  to  be  independent  of  speed,  in  the  sense 
that  the  same  benefit  is  received  for  a  particular  rudder 
sequence  despite  the  initial  speed.   For  example,  if  a  20 
degree  rudder  sequence  reduces  snap  roll  by  1^%   for  u   =20  kts, 

then  the  same  sequence  will  also  reduce  snap  roll  by  15%   for 

a  u  =  30  kts. 
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III. k      SAIL  SIZE 

A  submarine  sail  produces  a  large  roll  moment  which  is 
caused  by  the  sideslip  velocity.   The  magnitude  of  the  roll 
moment  varies  widely  from  submarine  to  submarine,  primarily 
due  to  the  extensive  variation  in  sail  location  and  size. 
The  sail  dimensions  are  primarily  derived  from  considerations 
for  the  quantity  of  equipment  and  the  dimensions  of  the 
equipment  to  be  placed  in  the  sail.   With  the  exception  of 
the  fairwater,  which  encloses  the  sail,  the  naval  architect 
pays  little  attention  to  the  hydrodynamic  considerations  of 
sail  design.   It  was  felt  that  if  the  sail  size  could  be 
reduced,  that  a  corresponding  reduction  in  snap  roll  would/ 
be  realized. 

An  investigation  was  performed  to  determine  if  the 
chord  or  the  span  had  the  greater  impact  on  snap  roll. 
Figure  J>-6   shows  the  results.   Obviously,  a  change  in  the 
span  will  have  a  greater  impact  on  roll,  as  evidenced  by 
the  larger  slopes  for  the  three  coefficients  considered. 
This  result  was  anticipated  since  the  moment  arm  associated 
with  the  spanwise  location  of  the  center  of  pressure  will 
have  a  greater  influence  on  roll  than  the  moment  arm 
associated  with  the  chordwise  location.   Additionally,  a 
reduction  in  span  will  produce  a  corresponding  reduction  in 
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the  lift  due  to  the  lower  aspect  ratio.   A  decrease  in  the 
chord  would  obviously  increase  the  aspect  ratio.   Conse- 
quently, the  span  was  reduced  by  1,  5»  10,  15,  20  and  30 
percent  for  a  number  of  submarines.   The  results  from  these 
computer  runs  will  be  presented;  however,  first  the  method 
used  to  predict  the  impact  on  the  hydrodynamic  coefficients 
of  reducing  the  span  will  be  presented. 

The  approximate  method  for  determining  the  impact  of 
the  sail  on  the  hydrodynamic  coefficients  was  partially 
derived  from  Abkowitz's  presentation  for  a  lifting  foil  in 
reference  (4-).   This  approximate  method  involves  determining 
the  Y  force,  N  moment  and  K  moment  at  the  sail  as  a  result  of 
velocity  disturbances  v,  r,  and  p.   These  coefficients  are 
reduced  correspondingly  for  a  reduction  in  the  sail  size. 
The  new  sail  coefficients  are  added  to  the  original  co- 
efficients, less  the  original  sail  contribution.   Then  the 
equations  of  motion  are  solved  for  r.      Figure  3-7  defines 
the  coordinate  system,  the  variables  involved  in  the 
derivation,  and  sign  conventions. 

For  most  purposes,  the  sail  x  and  y  velocity  components 
are  the  same  as  the  submarines' s  x  and  y  velocity  components, 
i.e. 


u  -  u      v  =  v 
s  s 
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The  subscript  "s"  denotes  the  local  condition  at  the 
sail.   The  angle  of  attack  produced  at  the  sail  by  the 
transverse  velocity  is 


ft  -   tan"1  Is  =  t 


-1  v  _  j.   -1  v 
^an 


Consequently,  we  can  calculate  the  Y  force  and  N  moment 
produced  at  the  sail  as 

Y  =  -  (L  cosy3  +  D  siryS) 


N   =  Y   .  x_ 

s    s    s 


where  L  and  D  are  the  lift  and  drag  forces  on  the  sail, 
respectively.   In  terms  of  their  respective  non-dimensional 
coefficients,  the  lift  and  drag  are 


L  *  GL  •  i   AsUs2  =  CL  •  Wus2  +  Vs2) 


D  »  CD  .  i/As(us2  ♦  vs2) 


where  A  is  the  projected  sail  area  and  e   is  the  fluid 
density. 

From  airfoil  theory  (6),  the  slope  of  the  lift  co- 
efficient curve  versus  the  angle  of  attack  is  shown  to  be 
approximated  by 
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AR  is  the  "reflected  aspect  ratio."   This  is  equal  to  twice 
the  mean  span  divided  by  the  mean  chord,  which  corrects  for 
the  proper  pressure  distribution  on  the  sail.   Using  the 
above  relations,  it  can  be  shown  that  the  Y  force  on  the 
sail  due  to  a  lateral  velocity,  v,  is 


(Y^  =  -*V-VH&n  *  V 


and  the  N  moment  due  to  v  is 


<Vs  "  <Vs  •  xs 


Similarly,  for  the  K  moment 


<Kv's  "  <Vs  •  *s 


For  the  effect  of  an  angular  velocity,  v,  where 

v  =  x   .  r,  we  find 
s    s 


Ys  due  to  r  =  (Yy)s  .  *„  .  r 

2 
N  due  to  r  =  (Y„)   .  x„   .  r     and 
s  v  s    s 


Ks  due  to  r  =  (Yy)g  .  xg  .  zs  .  r 
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Taking  the  derivative  of  these  coefficients  with  respect  to 
r, 


(Y  )   =  (Y  )   •  x 
v  r  s     v  s    : 


(N  )   =  (Y  )   -x 
v  r  s     v's    s 


2 


(K  )   =  (Y  )   -x   •  z 
r  s     v  s    s    : 


A  similar  process  for  the  contribution  of  an  angular  velocity 

with  v  =  -z   •  o,  yields 
s     s   *'   J 


(Y  )   =  -(Y  )   ■  z 

v  p's      v's    s 


<Vs  "  "<Vs 


<Vs  =  -<V 


A  similar  derivation  provides  the  hydrodynamic  coefficients 
for  the  sail  with  respect  to  the  accelerations  v,  r  and  p. 
The  key  coefficient  is,  of  course,  (Y- )  .   Based  on  the  added 

V   o 

mass  calculation  for  a  flat  plate  of  dimensions  s/2  and  c, 
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This  form  of  (Y- )   is  used  with  a  correction  factor  de- 
veloped by  Abkowitz  (197?)  and  finally  becomes 


(Y-) 
x  v  s 


"*  ~  z 


Kfs2  c2 


rp 


^s2  +  cZ 


1  - 


0.  5^ 


(1  +  s/c  +  c/s) 


where  s  is  the  mean  "reflected  span"  and  c  is  the  mean 

chord.   The  mean  "reflected  span"  is  directly  related  to 

the  "reflected  aspect  ratio",  in  that,  it  corrects  for  the 

proper  pressure  distribution  on  the  sail. 

Since  x  and  z   are  the  distances  from  the  center  of 
s      s 

pressure  to  the  center  of  gravity  along  their  respective 
axes,  the  following  relations  (7)  proved  beneficial, 


(CP)-  =  0.25  •  c 


(CP)i  =  (V3)7T-  s 


(CP)-  is  the  distance  to  the  center  of  pressure  chordwise 
location  from  the  sail's  leading  edge.   Similarly,  (CP)- 
is  the  distance  from  the  root  chord  to  the  center  of 
pressure. 

Employing  the  above  described  procedure  and  equations, 
new  coefficients  were  developed  for  reduction  in  sail  size 
for  a  number  of  submarines.   The  results  presented  in 
figure  3-8  are  for  two  different  submarines,  with  what  are 
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termed  large  and  small  spans.   "Large"  and  "small"  relate  to 
the  relative  size  of  the  sails  when  compared  to  similar 
submarine  designs.   It  can  be  seen  that  a  greater  benefit 
is  realized  by  the  larger  sail.   The  maximum  reduction 
applied  to  the  span  was  restricted  to  30fo.      Any  larger 
reduction  in  span  would  not  be  feasible  due  to  physical 
restrictions,  such  as  minimum  periscope  height  requirements. 

The  decrease  in  the  maximum  roll  angle  through  sail 
size  reduction  was  not  as  large  as  was  expected.   It  is 
believed  that  this  may  be  explained  by  the  fact  that  the 
sail  contributes  to  the  roll  damping  at  the  same  time  that 
it  is  acting  as  a  roll  producer.   Therefore,  by  reducing 
the  sail  size,  the  roll  damping  is  reduced,  as  well  as  the 
roll  production.   Fortunately,  as  indicated  by  figure  3-8 1 
roll  damping  appears  to  decrease  less  rapidly  than  roll 
production  for  a  given  sail  size  reduction. 

III. 5   SPEED  REDUCTION 

It  has  been  well  established  that  the  forward  velocity, 
u,  has  a  strong  influence  on  the  roll  angle  when  in  a  turn. 
It  can  be  shown  that  the  maximum  roll  angle  is  proportional 
to  the  square  of  the  velocity.   Therefore,  it  is  expected 
that  small  speed  reductions  will  give  significant  reductions 
in  ^   .   This  is  an  important  point,  since  if  a  submarine 


^2 


is  designed  for  a  maximum  speed  in  excess  of  30  kts ,  then 
the  operating  community  would  like  to  be  able  to  take  full 
advantage  of  this  speed  in  various  maneuvers ,  as  well  as 
in  a  straight  line  motion.   The  greatest  difficulty  with 
snap  roll  occurs  at  a  submarine's  top  speed.   It  would  be 
attractive,  from  an  operational  point  of  view,  to  be  able 
to  enter  a  turn  at  almost  maximum  speed.   With  this  in 
mind,  several  submarines  were  simulated  to  enter  turns  at 
from  one  to  four  knots  below  their  maximum  speed.   The 
results  are  presented  in  figure  3-9-   As  was  expected, 
small  changes  in  the  velocity  produce  relatively  large 
reductions  in  snap  roll.   For  example,  a  10$  reduction  in 
u,  decreases  the  maximum  roll  angle  by  almost  2,0%.      This 
method  of  reducing  the  effect  of  snap  roll  is  equally 
feasible  for  use  in  an  automatic  control  system  or  as  a 
manual  method  of  operation. 
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CHAPTER  IV  -  CONCLUSIONS  AND  RECOMMENDATIONS 

IV.  1   CONCLUSIONS  AMD  RECOMMENDATIONS 

From  this  study,  it  must  "be  concluded  that  we  can,  with 
sufficient  engineering  accuracy,  predict  full-scale  sub- 
marine dynamic  performance,  at  least  in  roll,  from  model 
data.   This  is  accomplished  through  the  use  of  a  computer 
simulation  model.   This  is  of  importance  to  the  submarine 
designer  for  three  reasons.   First,  this  prediction  tech- 
nique provides  the  ship  designer  with  an  invaluable  design 
tool.   The  stability  characteristics  of  a  multi-million 
dollar  submarine  design  can  he  accurately  predicted  prior  to 
making  a  major  commitment  to  the  procurement  effort.   The 
ship  designer  now  can  proceed  with  confidence  on  a  design, 
at  least  from  a  dynamic  stability  standpoint.   Secondly,  the 
computer  model  can  be  used  by  the  submarine  designer  to 
perturbate  and  iterate  on  a  particular  baseline  design. 
Consequently,  such  a  technique  could  be  employed  to  optimize 
a  design  with  respect  to  its  stability  criteria.   Finally, 
this  method  of  predicting  s.tability  performance  can  be 
employed  by  the  designer  to  establish  design  criteria. 
Without  a  foreknowledge  of  a  submarine's  performance  in  a 
particular  maneuver,  the  ship  designer  can  test  past 
similar  designs  to  establish  some  pattern  of  performance 
and  thereby  derive  a  stability  performance  criteria. 


It  is  apparent  from  the  results  of  Chapter  III  that  a 
small  additional  amount  of  GM  in  any  submarine  design  will 
significantly  reduce  snap  roll.   Bearing  this  in  mind,  it 
would  seem  appropriate  to  suggest  that  the  U.S.  Navy  revise 
its  present  design  criteria  for  GM.   The  present  criteria 
requires  a  minimum  metacentric  height,  which  does  not  work 
well  for  all  designs  from  a  stability  standpoint.   It  is 
proposed  that  the  revised  design  criteria  for  GM  be  based  on 
a  maximum  acceptable  snap  roll  angle.   This  would  be  a  more 
flexible  design  criteria  which  could  be  tailored  to  indi- 
vidual designs.   More  importantly,  it  would  couple  a  per- 
formance measure  to  the  naval  architectural  characteristics. 
It  would  also  accomodate  for  the  variations  in  size  of  the 
roll  moment  producers  which  are  encountered  from  design  to 
design,  i.e.  the  larger  the  roll  moment  producers,  the 
larger  the  required  GM  to  minimize  the  effect  of  snap  roll. 

The  results  of  the  other  methods  considered  in  Chapter 
III  suggest  obvious  measures  to  be  taken  in  an  effort  to 
reduce  snap  roll;  these  are  to  reduce  speed,  reduce  initial 
rudder  angle  and  reduce  sail  size.   As  was  indicated  in 
Section  III. k,    the  reduction  of  sail  size  was  not  as 
effective  in  reducing  snap  roll  as  was  anticipated.   The 
possible  reason  for  this  was  also  suggested  in  Section  III. 4 
However,  a  more  thorough  study  in  this  area  is  certainly 
warranted.   This  is  particularly  true,  in  view  of  the  fact 
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that  an  approximate  method  was  used  for  predicting  the 
appropriate  hydrodynamic  coefficients  for  the  contribution 
of  the  sail  to  snap  roll.   A  more  precise  method  of  pre- 
diction, coupled  with  a  careful  consideration  of  the  total 
effect  of  the  sail  and  sailplanes  on  roll  stability,  could 
lead  to  a  better  insight  into  the  causes  and  minimization 
of  snap  roll.   Other  recommendations  for  future  investi- 
gations are: 

•  The  use  of  differential  sailplanes  and/or 
sternplanes  as  roll  control  fins. 

•  The  use  of  unbalanced  rudders. 

•  The  effect  of  increasing  the  automation  of 
submarine  control  systems  on  stability 
characteristics . 

Three  areas  of  improvement  for  the  computer  simulation 
model  are  proposed.   The  first  is  the  extension  of  the 
present  model  to  use  the  six  degrees  of  freedom  equations  of 
motion.   This  should  not  only  improve  the  accuracy  of  the 
results  obtained  from  the  present  four  degree  of  freedom 
program,  but  it  v/ould  also  allow  for  the  investigation  of 
more  complex  maneuvers  and  dynamic  reactions,  e.g.  squat 
with  trim  plus  roll.   Secondly,  the  model  should  be  modified 
to  permit  for  more  ease  in  its  use  by  submarine  designers. 
The  less  complex  the  input  requirements  and  the  more  lucid 
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the  output,  the  wider  use  the  computer  simulation  model  will 
experience  within  the  design  community.   Lastly,  it  is 
recommended  that  the  present  program  and  future  versions  be 
optimized  with  respect  to  computer  time.   With  the  in- 
creasing expense  and  the  time-sharing  requirements  of  com- 
puter facilities,  wider  use  of  this  type  of  computer 
simulation  model  will  be  realized  by  an  efficient  use  of 
available  computer  time. 
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A.  2      COMPUTER    PROGRAM 
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A . k     INPUT/OUTPUT  variables 


VARIABLE 


MEANING 


FORMAT  DIMENSIONS 


INPUT/ 
OUTPUT 


AI,  BI,  CI 


Set  of  con- 
stants repre- 
senting the 
propellar 
thrust  in  the 
X-equation 


D10.4 


ND 


Al,  A2 


Limits  used 
for  selecting 
the  proper 
propeller 
thrust 


D10.4 


ND 


B 


Ship's 
buoyancy 


D10.4 


lbs 


COURSE 


Amount  of       DlO.'l 
ordered  course 
change 


deg 


I.  0 


DELB 


Initial  sail- 
plane 
deflection 


DIO.^J- 


deg 


DELGI" 


Amount  the 
original  GM 
is  to  be 
changed 


D10.4- 


ft 


DELRO 


Rudder  de- 
flection at 
time  T 


D9-3 


deg 


ND  =  non-dimensional 
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Initial 

sternplanc 

deflection 


DlO.'i    de; 


t> 


D2LT 


Time  increment  D10. k 

used  in 

iteration 


sec 


I,  0 


ID 


^-0  character 
alpha- 
numeric 
heading 


Al 


ND 


I,  0 


INDEX 


If  >0,  read 
new  submarine 
coefficients . 
If  ^0,  run 
same  sub- 
marine for  new 
initial 
conditions. 


12 


ND 


IX,  IY,  IZ 
IXY,  IXZ,  IYZ 


Moments  of 
inertia 


D10.^ 


ND 


KCOEFF 


K-equation 
coefficients 


D10.4 


ND 


Ship's  over- 
all length 


D10.4 


ft 


M 


NCOEFF 


Ship's  mass     DIO.^J-     ND 


N-equation      D10.^     ND 
coefficients 


I 
I 


71 


P,  Q,  R         Initial         DlO.^i   deg/sec 

angular 
velocity 
about  the  x, 
y  and  z  axes, 
respectively 


PDT,  QDT,  RDT   Initial         D10.4   deg/sec2      I 

angular 
accelerations 
about  the  x, 
y  and  z  axes, 
respectively 

PDTO,  RDTO      Angular         D9 • 3    deg/sec2      0 

accelerations 
about  the  x 
and  z  axes, 
respectively 


PHI  Initial  angle   D10./+     deg 

of  roll 


PHIO  Roll  angle  at   D9-3      deg         0 

time  T 


PSI  Initial  angle   D10.4     deg 

of  yaw 


PSIO  Yaw  angle  at    D9-3      deg         0 

time  T 


RHO 


Sea  water       D10.4-   slugs/ft^ 
density 


RO  Angular         D9 . 3    deg/sec       0 

velocity  about 
the  z  axis  at 
time  T 
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RRATE 


Average 
rudder  rate 


DIO.'J-        deg/sec  I,    0 


RUDAI'T 


Maximum 
ordered 
rudder 
deflection 


DIO.^    de.-r 


1,0 


RUD1 


Value  of 
reduced 
rudder 
deflection 
associated 
with  sub- 
routine ALTRUD 


D10./+ 


deg 


I,  0 


'resent  time    D9«3 


;ec 


THETA 


Initial  angle   D10.il-    deg 
of  pitch 


TLAG 


Time  lag  of 
of  the  control 

system 


D10.il- 


:ec 


U,  V,  W 


Initial 
forward , 
lateral  and 
vertical 
velocities 


D10.ii-    kts 


UDT,  VDT,  r'/DT 


Initial 
forward, 
lateral  and 
vertical 
accelerations 


D10.il    kts/hr 
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UDTO,  VDTO      Forward  and     D9 . 3    kts/hr        0 

lateral 
acceleration, 
respectively 
at  time  T 


IJO  Initial  D9-3  kts  0 

forward 


velocity 


U00,  VO         Forward  and     W-3  kts  0 

lateral 
velocities , 
respectively 
at  time  T 


WT  Ship's  weight   DIO.^    lbs  I 


XB,  YB,  ZB      The  x,  y,  z     D10.'+    ft  I 

position  of 
the  center  of 
buoyancy 


XCOEFF  X-equation      DIO.'J-    ND 

coefficients 


XG,  YG,  ZG      The  x,  y,  z     DIOJl-     ft 

position  of 
the  center  of 
gravity 


XX,  YY  The  moving  x    D9-3     yds  0 

and  y  coordi- 
nates of  the 
submarine' s 
origin  v/ith 
respect  to  a 
fixed  coordi- 
nate system, 
whose  origin 
is  the  point 
where  the  turn 
begins 

7^ 


YCOEFF  Y-equation      D10.4     MD 

coefficients 


15 


APPENDIX  B 


B.l   AXIAL  FORCE 


m 


2     2 
u  -  vr  +  wq  -  xr    (q   +  r  )  f  yp(pq  -  r)  +  z^(pr  +  q) 


'(] 


G 


/  /''  X   'q2  +  X   'r2  +  X   '  rp] 
2  A  [_qq  H  rr       rp  n 


+  X   '  vr  +  X   'wq 
vr       wq 


2 


D2  X   'u 


2         2         2 
^  +  X   'v   +  X   'w^ 

vv       ww 


4i 


2  2 


2 


2  * 


HrSr^   +X<fsSs'^2   "  %  ft/& 


2  2 

a.u'    +   b.uu      +   cu 
1  1       c  1    c 


-    (W   -  B)    sin  9 


.r2u2 

r 
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B.2   LATERAL  FORCE 


m 


2  2 

v  -  wp   +  ur   -  yG(r      +   p")    +   zQ(qr   -  p)    +  xQ(qp   +  r 


r)J    = 


'hi 


k 


Y-'r    +   Y-'p   +   Y     .    :  p  I  p  J     +   Y      'pq   +•  Y      '  qr 

r  P  P  |P|         '  Pq  qr  H 


+  li 


3 


Y- 'v   +   Y      ' vq    +   Y      'wp  +   Y 
v  vq  wp      1  wr 


'  wr 


fli 


Y    'ur   +  Y    'up   +  Y,    ,c    'ulr  |£ 
r  P  |r|or      '     '  ui 


+   Y    ,    ,  'v  |(v2    -  w2)~2~  I  |r 
v'rl   Ivl 


*li 


2  2,1 


Y*'u      +   Yv'uv   ♦  Yv|v|'vl(v      +  w    )2' 


+  £ 


2 


JT^vw'™   +  Y<Jr'u2<5r] 


+    (VI    -   li)    cos    8   sin  f 


*  fjf  Yr?'ur(7-    1) 


+  !l 


2  2n4 


Y      'uv   +   Y     .    Im  'v  |(v2    +  w2)2  I    +   Yf      *J 


.u 


(7  -   1) 
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B.3      ROLLING  MOMENT 


V   +    (Iz    "    V^    "    <**    +    ^)Jxz    +    (r2    "   ^2)Iyz    +    (^r    "   ^xy 


+   m 


yQ(w    -   uq    +   vp)    -    zG(v   -   wp    I   ur) 


2  * 


K-'p   +   K-'r    t   K      'qr    +   K      ' pq    +   K     ,    ,'olp 


qr 


pq   F4    T      P|p| 


+  f* 


/|. 


K    'ud   f  K    'ur   +  K- 
p     *    '      r  v 


.,i 


*ii 


4 


K      ' vq    t   K      'wp   +  K 
vq      H  wp      ^  wr 


'  wr 


fliT 


*'u2    +  Kv'uv   +  Kv.|'v|(v2   +  w2)2 


•i/! 


K      'vw    •»■   Kr    ' 
vw  Or 


u2£rl 


1-    (yrW   -  yRB)    cos    9   cosf  -    (zpW   -   zRB)    cos    9   sin  f 


G"         JB 


G"  BJ 


3^       .,?- 


+  /^K*fu  (?-  i) 
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B.4      YAWING  MOMENT 


V   +    (Iy   "   V***   "    (*   +   rp)Iyz    +    ^   -   P2)Ixy   +    ^   '   -^!zx 


+   m 


xG(v   -  wp   +  ur)    -  yG(u   -  vr   +  wq) 


li 


N-'r   +  N-'p  +  N      'pq   +  N      'qr   +  N,    i  '  r  |r  | 


pq 


qr  H  r  Jr| 


+  li 


Ly 


N  •  '  v   +   N      '  wr    \   N      '  wp   +   N      f  vq 
v  wr  wp      ±  vq      ^ 


♦|i 


4 


N    '  up   +  N    'ur   +  N,    |  r    'u|r|ir   +   N,    ,     '  I 
p     l  r  |r|dr  '  |v|r    ' 


'  |(v2   4  w2) 
v  r 


*|r] 


*iJ 


.,,2 


9     1 


N*'u^    t  Nv'uv   +   N-.'v  Kv      +  w^)2| 


♦sri 


N      *vw   +  Nr   '  u25r 
vw  or 


] 


+    (xr  W   -   xp   B)    cos    0   sinf+    (yp  W    -  yR  B)    sin   9 


'G 


G  "        'B 


4  |^Nr?'ur    (7[ 


-   1) 


>u 


3 


N       'liv   +   N     ,     i-'v  |(v2    +   w2)2  |    +   Nr      *  X  u2 


(7  -  1) 


79 


Thesis  17' 2  63 

1.625  Libera  to  re 

An    investigation   of 
snap   roll    in   submarines. 


4     OCT  8d 


7  5 
9  7  977 


Thes  is 
L625 


171169 
Li  beratore 

An    investigation  of 
snap    roll    in   submarines. 


thesL625 

An  investigation  of  snap  roll  in  submari 


1 

tjlj 

ijMIIH 

II 

1 

3  2768  002  11759  0 

DUDLEY  KNOX  LIBRARY 


